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SUMMARY 

Evaluation and analysis were made of the mean and turbulent terms 
of the equations of motion and the stress tensor at four stations in a 
turbulent boimdary layer with a progressively increasing adverse pres- 
sure gradient. 

Good agreement between the values of skin friction obtained by heat 
transfer - skin friction instrument measurements and by the evaluation 
of the Ludwieg-Tillmann empirical equation was found to exist. The eval- 
uation of skin friction from the integrated momentum equation failed to 
agree with the results of other methods of obtaining skin friction. 

Evaluation of the terms of the tixrbulent stress tensor indicates 

that the normal stresses pv^ and pw^ and the shear stress -puv are 

of comparable magnitude, while the normal stress pu^ was found to be 
roughly four times as large near the wall. The angle between the prin- 
cipal £LXis of the turbulent stress tensor and the boundary appeared to 
be largely independent of x-distance or pressure gradient . 

Evaluation of the x- and y-direction equations of motion shows that 
the rate of change normal to the wall of the mean square of the y- 
direction turbulent velocity in the y-direction equation when taken near 
the wall is as large as or larger than any term in the x-direction equa- 
tion. The X- and y-momentum equations are, however, still independent. 


IlttRODUCTION 

In recent years progress has been made in the measurement of the 
statistical properties of specialized turbulent shear flows. For the 
free turbulent shear flows in the wakes of circular cylinders and in 
free air jets, Townsend (ref. l) and Corrsin (ref. 2), respectively, 
have evaluated from experiments the important statistical terms. The 
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most accurate investigations have been made in free turbulent shear flows, 
g absence of solid boundaries greatly facilitates measurements* 

In fully developed t\irbulent flow in pipes and channels, which is the 
other extreme, Laufer (refs. 3 and 4) has presented extensive measure- 
ments. Laufer's results brought to light significant features of flows 
with solid boimdaries; it was foimd that flow conditions in the close 
proximity of the wall were of major importance. For turbulent boimdary 
layer flows, which represent an intermediate type of flow between fully 
free turbulent and channel flow, measurements are meager. The very 
recent work of Townsend (ref. 5) and Klebanoff at the National Bureau of 
Standards has made inroads into the understanding of the characterist ics 
and structure of turbiilent boundary layers with zero pressure gradient, 
with both investigators directing their attention toward the turbulent 
energy balance. The exacting investigations of Askenas, Riddell, and 
Rott at Cornell University in a zero pressure gradient boxmdary layer 
were imable to eliminate the apparent inconsistencies in the measured 
shear distributions. (The shear stress curves possessed a maximum away 
from the wall.) Preliminary experimental measurements with hot-wire 
equipment by Schubauer and Klebanoff (ref. 6) in boundary layers with 
adverse pressure gradients also resulted in inconsistencies in the shear 
distribution. (The shear stress curves indicated too high a value of 
skin friction.) Other experimenters have also been confronted with in- 
consistent wall shear stress evaluations. 

The NACA Lewis laboratory is employing a 6- by 60-inch boundary 
layer channel in a long-range program to provide an insight into the 
phenomenon of turbulent boundary layer separation. The method of attack 
is to evaluate all mean and fluctuating quantities of interest (a) far 
upstream in the unseparated turbulent boundary layer, and (b) within 
the region of separation. A comparison of the two sets of data should 
then greatly assist in the study of the manner in which the boundary 
layer separates . 

This report constitutes the first phase of part (a), in which all 
terms pertaining to momentum have been evaluated, such as the mean and 
turbulent shear stress and stress tensor, along with the time-averaged 
equations of motion. These terms are evaluated and analyzed for a flow 
with progressively increasing pressijre gradient. Special effort was 
made to present a consistent set of measurements j however, wherever con- 
fliciting results were encountered, two or more independent sets of data 
are shown. 


THEORETICAL BACKGROUND 

poliowing the assumptions of Reynolds, the txme— averaged Navier— 
Stokes equations of motion for an inconipressible turbulent fluid motion 
may be written 


3173 


CF-1 back 


MCA TW 3264 


3 



( 1 ) 


where 


(i,j = 1,2,3) 
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and 
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(a list of symbols will be found in appendix; A.) The two basic groups 
of terms appearing in equation (l) are the inertia terms, which are 
those on the left side, and the stress terms, which appear in the pax- 
entheses on the right side. 

For the particular boundary layer flow studied in this report, 
equation (l) may be simplified by assuming that the mean motion of the 
flow is two-dimensional; thus all gradients in the z-direction are zero 
Although the tiirbulent motions are of a three-dimensional natiare, both 
they and the mean flow are assumed symmetric relative to the xy -plane. 
From this condition of plane -symmetry, it follows that a reflection of 
the z-axis will require that uw and vw equal their respective neg- 
atives, which implies that both must vanish. This argument for the 
vanishing of uw and w is patterned after the argument employed in 
isotropic turbulence, where the turbulence is spherically symmetric 
(ref. 7). Consequently, for the flow investigated, the steady-state 
equations of motion may be simplified to, for the x-direction. 
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and for the y-direction 
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Previous measiirements (ref. 6) of the terms in equations (2) and 
(3) have not been sufficiently accurate to discuss the equations ade- 
quately, but it is generally believed that certain of the terms may be 
neglected, thus further simplifying the equations. 


The stress terms for plane -symmetry tixrbulent shear flow may be 
arranged as a symmetrical stress tensor . 
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These terms are the flux of momentum in the flow and their deriva- 
tives yield the forces acting within the boundary layer. It will be con- 
venient to further subdivide the stress tensor into two additive tensors, 
one composed of the stresses associated with viscosity and pressxire and 
the other containing the turbilLent contributions. 

The specific object of this report was the experimental evaluation 
of the stress terms, the eqioations of motion, and other terms that are 
known to be important for the understanding of turbulent boundary layer 

flow. 


TEST PROCEDURE AND RESULTS 

Tunnel and test conditions . - The measurements reported herein were 
made in the turbulent boundary layer along the test wall of the Lewis 
6- by 60-inch boundary layer channel. A schematic diagram of the tunnel 
is shown in figure 1, and a complete description of the channel, except 
for a modified inlet, appears in reference 8. All measurements were 
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taken -with a constant Reynolds number per foot of 3.33x10^^ which cor- 
responded to a free -stream velocity of approximately 58 feet per second^ 
maintained at the entrance of the test section. The suction across the 
porous wall was sufficient to keep the boundary layer along the porous 
wall at a constant thickness (the same as reported in ref. 8 for the 25 
in. water pressure drop case). 

The four stations at which the momentum surveys were tak;en are shown 
in figure 1. They are approximately 18 inches apart, with the first 
located 72 inches from the inlet contraction. The boundary layer was 
approximately 1 inch thick at station 1 and increased to 3 inches at 
station 4. 

Pressure measurements . - The static -pressure distribution along the 
test wall was obtained from wall static orifices (0. 025-in. -diam. ) loca- 
ted on lines 12 inches above and below the center line, as noted in fig- 
ure 2. The actual pressures were recorded with a water micromanometer 
with a least count of 0.001 inch of water. The manometer was referenced 
to atmospheric pressure. 

Mean velocity profiles were evaluated from measirrements with the 
total -pressure probe shown in figure 3. The static pressure was obtained 
from wall orifices. The probe tip was of 0.040-inch-outside-diameter 
tubing with 0 .004 -inch-thick walls. The end was flattened to give a 
rectangular opening with an over-all height of 0.020 inch, thus allowing 
the center of the probe opening to be brought within 0.010 inch of the 
wall. The pressures were recorded using a ±0.15 psi and a ±0.05 psi 
full range strain-gage -type pressure transducer. The transducers are 
of the differential -pressure type consisting of a bellows and an unbonded 
strain gage. The output of the transducers was recorded on a 2 millivolt 
full-scale recording potentiometer, so that under steady conditions it 
was possible to detect pressure differences of the order of 0.01 and 
0.004 pound per square foot, respectively, for the two instruments. The 
voltage balance and calibration circuits used to operate the transducers 
are shown in figure 4. 

For each profile the probe was zeroed at the wall by an electrical 
contact and then actuated outward through the boundary layer with an 
actuator which could be set at any distance in the boundary layer with 
an error no greater than ±0.001 inch. 

It is of interest to point out a physical phenomenon which has been 
encountered in all turbulent boundary layer measurements made with both 
pressure transducers and manometers. Figure 5 shows inserts of the out- 
put signal of a pressure transducer connected to the total -pres sure 
probe at different y-distances within the boundary layer of station 1. 
These records show large -amplitude fluctuations occurring near the bend 
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in the velocity profile; with greatly damped signals being noted near 
the wall. This fluctuation record should not be inteipireted quantita- 
tively, since the output of the system will be influenced by the natural 
frequency (250 cps) of the pressure transducer and by the frequency re- 
sponce of the recording instmoments, but it does give a qualitative idea 
of the behavior to be expected when total -pres sure measurements of this 
type are made. The fluctuating signal was shunted out through a capa- 
citance in order that mean pressure could be recorded. It was noted 
that the indicated mean pressures were higher than the approximate mean 
of the fluctuating signal. 

The calculated mean velocity profiles are shown in figure 6. In 
all cases except the last station (station 4), the recording equipment 
could be sufficiently damped to obtain \miform velocity profiles. For 
station 4 it was necessary to obtain the points near the wall as the 
average of several points taJcen at different times for each y-distance 
because of the large fluctuations. Hot-wire velocity measurements were 
also included in the averaging of the points near the wall for station 4. 

The mean velocity evaluated from hot-wire measurements has also 
been included in figure 6(b). Although the hot wires employed for fluc- 
tuation measurements are rather poor for mean velocity measurements be- 
cause of calibration shifts caused by dust particles striking the wires, 
the selected results presented in figure 6(b) are in good agreement. 


The correction suggested in reference 9 for the effect of turbu- 
lence on the total -pres sure measurements has been evaluated and the cor 
rected mean velocities at a few points are included in figure 6. The 
effect due to t\irbulence results in a maximum error for the mean veloc- 
ity of roughly 2 percent which is for the region near the wall. The 
correction for the displacement of the effective center of the probe 
opening due to the velocity gradient across the mouth of the probe ^ref 
10 ) was also included in the corrected points. 


Flow deflection measurement. - The flow deflection angle measuring 
device is shown in fi^e 7. The instrument employs the standard two 
total-pressure tubes, each making an angle of 45° with the mean flow. 

It is a null -balance system, vith a special actuation mechanism being 
employed to rotate the probes in the boundary layer. The accuracy o 
measuring the flow deflection was of the order of lO.l - The ^gles 
measured between the mean flow and the test wall are shown in fi^e 8. 
The first station proved to have too small a variation in deflection 
angle in the boundary layer; thus no measurements could be obtained. 
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Heat transfer - skin friction measurements . - The technique of de- 
termining the local turbulent skin friction from a measure of the local 
heat transferred from a small locally heated insulated segment of the 
wall was developed in reference 11. For the present report a calibrated 
relation between the skin friction and the heat transfer was obtained in 
a fully developed turbulent channel flow^ where the skin friction coiild 
be independently determined from the pressure gradient. A description 
of the instruments employed and the calibration technique used is given 
in reference 8. 

Preliminary measurements (ref. 8) with the heat transfer - skin 
friction instruments gave only qualitative agreement with the Ludwieg- 
Tillmann predictions. However^ it was found in the present work that 
the agreement was improved by taking several measures of the heat trans- 
fer at a station and averaging the results to minimize random errors in- 
herent in the measure of such low voltages . The local measurements of 
skin friction obtained with the heat transfer - skin friction instrument 
are shown in figure 9. 

The simplicity of operation of these heat transfer - skin friction 
instruments makes them a very practical tool for the type of boundary 
layer studies now under way. 

Hot-wire anemometry measurements . - The hot-wire anemometer system 
employed for the measurement of turbulent fluctuations was the constant - 
temperature system described in detail in reference 12^ with minor changes 
made to reduce as much as possible the noise level of the modified system 
so that it now compares favorably with the conventional constant -current 
systems. The constant -temperature system reduced the time required to 
record data to less than half that for constant -current systems^ and di- 
rect comparisons show that the two systems give equally accurate measure- 
ments of the turbulence for the levels encountered. 

A detailed description of the system and the theory of measurement 
of fluctuating velocities by constant -temperature hot-wire anemometer 
systems is presented in appendixes B and C . 

The four particular types of hot-wire probes used for the measure- 
ments are shown in figure 3. The wires used were 0.0002-inch-diairieter 
etched tungsten (shc^xi as sample fig. 18^ ref. 12). The wires were 
mounted by copper plating the ends^ so that it is possible to soft- 
solder the wires to the probe prongs (see ref. 13). All wires had a 
0.040-inch unplated center section^ with the probe prongs being 0.050 
inch apart . The x-probes have their wires mounted in parallel planes 
0.010 inch apart ^ since this appealed to be the minimum distance for 
which practical mounting could be accomplished. The wire angles (45^ 
and 30^ with the direction of flow) were maintained within approximately 
±1^ with the use of a toolmaker's microscope. 
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In the operation of x-prohes no attempts were made to match hot 
wires or in any way to rotate the probes to match wire sensitivities, as 
there is some question as to the accuracy of using such methods in sim- 
plifying data recording. All measurements with x-prohes presented herein 
were made up of two or more independent readings of hot-wire outputs. 

The y-distance for the hot-wire probes was determined in the same 
manner as that described for the total -pressure probe, with extreme cau- 
tion being used to prevent more than a light contact with the wall. The 
construction of the probe limits the nearness to the wall for wh^h in- 
tensities may be measured to the following approximate values: u ^ 0.001 

inch, 0.010 inch, and v^, 0.025 inch. Also, the values of w^ near 
the wall are questionable, as the_two x-wires are in planes 0.010 inch 

apart; however, measurements of u^ with the w^ probes close to the 
wall check favorably with those obtained with a single-wire probe, so 
that the error is believed small. 


The turbulent intensity measurements for the four stations are pre- 
sented in figures 10 to 12. All measurements were reduced to the sta- 
tistical velocity fluctuations by using equations (,B5J, (.OllJ, and 
in appendixes B and C. Stations 1 and 2 show little random variation in 

measurements of while stations 3 and 4 indicate some random scat- 

ter- therefore, check measurements are included for stations 3 and 4 
(figs. 10(c) and (d) ) . A check of the consistency of the fluctuation 

measurements is included in figure 10(a), where the measurement of ^fu^ 
obtained by three independent hot -wire -probe types is shown. In all 

measurements of and present^, the maximum allowable ran- 
dom error in the calculated values of from the x-wires was re- 

quired to be less than 10 percent; any profiles indicating greater disa- 
greement were not included in the report . 


The ^ u^ measurements at station 3 are somewhat questionable, 
since considerable difficulty was encountered in obtaining co^istent 

measurements. A systematic disagreement was noted in the pro- 

files from day to day, with checks indicating that two distinctly dif- 
ferent types of profiles were equally likely. This duality of 

profiles will be discussed more fully in the section "Summary of velocity 
measurements . " 

The turbulent shear stress values obtained with a 45° x-probe are 
presented in figure 13. The meas\irements were reduced to stress values 
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^ by using equation (C6) in appendix C. In the course of these measure- 

ments it proved difficult to obtain consistent values of uv without a 
< great deal of scatter^ since the evalixation of uv requires working 

with the squares of the physical measurements of the hot-wire outputs, 
thus magnifying any error which might have been present in the actiial 
measurements. Figure 14 shows a set of points of the squared output of 
each wire of a uv probe; the faired curves from plots similar to this 
were used to obtain the faired cirrves shown in figure 13. The data 
points in figure 13 correspond to the direct measured difference between 
the hot wires and exhibit random errors as great as 30 or 40 percent. 

A consideration of the errors entering the measurement and calcu- 
lation of the turbulent velocity fluctuations along with the cross checks 
made on the data indicates the maximiom error for the measurements pre- 
sented (other than uv) is of the order of 10 percent. 


DISCUSSION OF RESULTS 

In order to obtain a picture of the momentum distribution through 
t\irbulent boundary layers, the terms of the modified Navier-Stokes equa- 
tion for turbulent flow and of the turbulent stress tensor will be eval- 
uated. However, to define the turbulent boundary layer flow investiga- 
ted, terms not necessarily associated directly with momentum, such as 
the mean flow parameters, the turbulent velocities, and the wall shear- 
ing stress, must also be evaluated and discussed. 

Mean flow parameters . - The mean flow parameters 6 , 0, Ren, and 
H evaluated from the faired mean velocity profiles axe presented in 
figures 15 and 16. Also included in the figures are points obtained 
when corrections for turbiilence (ref. 9) and probe displacement (ref. 

10 ) were applied to the data. The corrected values of 6 also include 
momentiim contributions due to the turbulent motion as derived ih refer- 
ence 14. Actually, the largest corrections to 0 were due to probe 
corrections and only secondary corrections were realized from the tur- 
bulence. In all, it appears that no great error is incurred in 0 by 
neglecting all corrections. The displacement thickness 8* shows a 
slightly greater effect due to probe corrections, with a similar effect 
on the profile fomn parameter H. 


Summary of velocity measurements . - Faired curv es o f the m ean and 

turbulent velocities are shown in figure 17. The and com- 

ponents of velocity were found to be very nearly equal over most of the 
boundary layer. The difference exceeds the probable error of measure- 
ment only near the boimdary. Except near the wall, f2 

and 17 

are seen to be roughly one -half 1^. 
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Some interesting points are noted in comparing the shapes of t he 

ttirbuient velocity profiles. The boundary conditions imposed on *\/~v2 
by the wall affect the distribution of v ^ relatively long distance 

out into the boundary layer, since the curves begin to decrease 

before the rate of energy converted to turbulence from the mean flow has 
reached a maximum^ as will be shown later (fig. 29 ) . Th is decrease in 

the ^ curve appears to affect the if^ and curves in that 

inflection points in both curves appear to occur at about the same point. 

It is believed that the conditions imposed on by the wall may not 

require the fluctuating energy to be dissipated directly, but may rather 

result in a reorientation of enerp- in the other two^rections . 

This concept of the shifting of the energy to the y and y 

components can be interpreted crudely as requiring that the axis of ro- 
tation of the turbulent eddies be turned more toward the perpendicular 
to the wall, as might be expected from the model propose^in reference 

15. From this concept it appears that near the wall y of the three 
turbulent components, is m ost directl y re lated to the boundary cona- 
tions of the flow, with and in turn depending on \v^. 


The variation of the profile with boundary layer development 

and pressure gradient also appears of interest . The general shapes of 

the \fu^ curves for the first two stations (figs. 17(a) and (b)) are 
very similar to those observed by Laufer in a f\illy developed turbulent 
charnel flow (ref. 3). The last station (fig. 17(d)), however, shows a 
decided change in the shape of the curve, with a corresp^ding change in 

the curve. In figure 17(c), the two possible profiles 

observed at station 3 for different runs are again shown. Of the two 
profiles, the solid line profile (l) agrees in form with profiles of 
the first two stations and the remaining profile agrees in form with 
that of the last station. 

A comparison of the mean velocity profiles shows much the same re- 

s\ilt as the v u^ profiles, in that the first three stations show a 
rather universally developing profile and the last station shows the 
development toward an inflection in the profile. While the actml phys- 
ical change in velocity profiles between stations 1 and 3 is slight, 
there is a great deal of change between the profiles at stations 3 and 4, 
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The previous observations can now be employed to attempt an explan- 
ation of the duality of measurements observed for station 3. In general^ 

the observations of ^ u^ at station 3 fall into two distinct groups^ 
one set agrees qualitatively with the profiles for the first two stations^ 
and the second agrees with the observations for the last station. It 
therefore appears that the flow conditions were not completely reproduced 
from day to day, although no unexplained inconsistencies were noted for 
the other three stations. The fluctuation in setting up flow conditions 
could not have been large, as inlet conditions were nearly the same from 
day to day, with corrections being made for changes in atmospheric pres- 
sure and temperature (analysis shows the free-stream velocity was changed 
by a maximum of 0.7 percent as inlet conditions were varied from day to 
day). This has led to the conclusion that the change of velocity dis- 
tribution to a profile as observed at station 4 must occur over a very 
short distance and is quite sensitive to flow conditions. The impression 
is that the profile changes abruptly from one type of distribution to^ 
another. Reference 16 contains some interesting information on equilib- 
rium profiles in turbulent boundary layers, which appears to indicate 
similar results. 


Evaluation of wall shearing stress . - Several methods are now avail- 
able for the evaluation of the wall shearing stress; hence it was desir- 
able to study the methods and atten 5 )t to place the necessary limitations 
on each. Figure 9 is a summary of the results of the several methods 
employed. The empirical equation of Ludwieg and Tillmann (ref. 17) 


= 0.246x10 


-0.678H -0.268 


R 


e 


and the empirical flat-plate equation of Falkner (ref. 18) 


( 5 ) 


, _ 0.01306 

(Reg)V6 


are compared with the values obtained from the heat -transfer measure- 
ments. As was foimd in the previous work (ref. 8), the heat -transfer 
measTirements and the Ludwieg-Tillmann relation are in good agreement. 
Secondly, the Falkner equation gives very good agreement for the first 
two stations and an error of less then 15 percent for station 3. The 
agreement of Falkner 's equation for the first two stations might have 
been expected, since figure 16 shows H to have remained nearly^ con- 
stant, while the last two stations show rather marked increases in H. 


The evaluation of skin friction by use of the von Karman integral 
momentum relation is also included in figure 9. The suggestions of each 
of references 14, 19, and 20 were considered in the evaluation of sk^ 
friction for the flows investigated in this report. Unfortunately, these 
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suggestions give negligil^le con'ti’ibu.tions over tlie values obtained, when 
a1 1 turbulence contributions are neglected. The resulting values of C^. 
obtained using the form suggested for 0.^. in reference 14 and using the 
wall pressure as suggested in reference 19 are shown in figure 9. 


It is evident that the suggested corrections have not been suffi- 
cient to allow the integrated momentum relation to predict skin friction 
in pressiire gradients, and it appears necessary to look elsewhere for 
the answer. Another possibility that has been suggested is the lack of 
two -dimensionality of the mean flow. One of the chief aims in selecting 
the dimensions of the 6- by 60-inch channel was to eliminate any end 
effects on the flow along the center. Velocity profile and static pres- 
sure surveys show no marked variation in the center region, so that the 
flow was assijmed to be two-dimensional in nature. However, the work of 
Clauser (ref. 16) in a tunnel at The Johns Hopkins University, which 
produces flow similar to that of the 6- by 60-inch channel, gave a strong 
indication of lateral spreading of the flow which resulted in an addi- 
tional loss of momentum not included in the two-dimensional momentum 
equation. 

It thus appears that for adverse pressure gradients the two- 
dimensional integral momentum relation is not adequate for evaluating 
skin friction. Moreover, most practical applications will involve even 
further deviation from two-dimensionality than those existing in the 
wind tunnel. 

The large discrepancies shown in figure 9 may also be due to the fact 
that the evaluation of skin friction from the momentum equation involves 
obtaining small differences of large quantities that are derivatives of 
experimental curves . 

Although lateral flow can cause large errors in skto friction cal- 
culated from the integral momentum equation, this flow is still so small 
that the actual skin friction and other physical features of the basic 
flow ane not appreciably affected. 


Shear st ress distribution across boundary layer . -To obtain the 
total stress distribution across the boundary layer it is necessary to 
combine the turbiilent and viscous contributions. Figure 18 shows the 
contribution of each of the shear stresses near the wall. Outside the 
region plotted, the viscous contribution will be negligible. 

Although it was impossible to measure the turbulent stresses closer 
than 0.025 inch from the wall, the trends indicate reasonable agreement 
with the Ludwieg-Tillmann values. Some indication of the distribution 
near the wall might be found, as suggested in reference 6, by evaluating 
the streamwise momentum equation at the wall. This yields the ini la 
slope of the total shear curve as 


_ 





( 7 ) 
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A line of this slope is included in figure 18 at the wall value predicted 
by the Ludwieg-Tillmann equation. 


There is also a limited amount of information available from figure 
18 regarding actual shear emd velocity distributions very close to the 
wall. The laminar sublayer has been assumed to be characterized by a 
linear velocity profile and vanishing turbulent shear stress for some 
small but finite distance out from the wall. Although actual data points 

do not extend deep into the sublayer^ the trends indicated by the p 

curves (fig. 18) suggest a parabolic rather than a linear velocity dis- 

constant). Equation (2) evaluated at the wall 


tribution (i.e.^ 




duv , 

- P 's — + 






h'l 


( 8 ) 


indicates that the assumption of a laminar sublayer with the conditions 
of a linear velocity profile could be consistent with the equation of 
motion only for a zero pressiire gradient. The assxamption often made that 
turbulent shear stress becomes iinimportant for values of yU^/v less than 
about 30 also appears inconsistent with the curves of figure 18, since 
the viscous shear stress terms indicate they reach the total shear stress 
veilue only very near the wall. It is noted that in the region next to 
the wall the total shear (viscous plus turbulent) changes very little 
compared with the large changes in the separate viscous and turbulent 
shear stress curves. It might be more usefiil, therefore, to speak of a 

dr 8p^ 

sublayer in which rather than of a "laminar" sublayer with 

linear velocity profile and zero turbulent shear stress. 

Attempts to predict the shear distribution through the turbulent 
boundary layer have been largely empirical in nature. Reference 21 
tries to reduce the prediction of shear distribution to a power series 
representation, but can obtain at most only a finite number of terms, 
since the evaluation depends entirely on limited boundary conditions. 
Reference 22 resorts to the boundary conditions and the equations of 
motion for a fully developed turbulent channel flow to predict a special 
power series for the velocity profile (which in turn is used to evaduate 
the shear profile). The power series of reference 22 results in good- 
agreement for the special one -dimensional flow investigated. Both 
methods must be classified as en 5 >irical, although reference 22 does rely 
on the equation of motion to give the final form of the shear distribution. 

A third method for predicting the shear distribution near the wall 
from mean velocity measurements is presented in reference 23. Although 
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this method must be limited to the vicinity of the wall, it appears to 
be of a general form applicable to any type of turbulent boundary layer 

flow. 


Starting with the assumption that the mean velocity profile near^ 
the wall can be represented by the universal profile predicted from sim' 
ilarity relations 


U_ 

Ut 



(9) 


which has been verified 
verified in this report 
are integrated. 


by many investigators in 
(fig. 19), the equations 


the past and is again 
of motion and continuity 


Equation (9) is used in the continuity equation and the bounda^ 
condition V = 0 at y = 0 is applied in reference 23 to obtain the 
following relation for V: 

V_ ^ iL (lO) 


The momentum equation was integrated, using equation (9), which with the 


boundary condition 
distribution 


T = T 


w 


at y = 0 gives an expression for the shear 


Z. 





( 11 ) 


where dp/dx has been assumed independent of y. These equations for 
T and Y depend only on the assumption of the existence of a relation 
as expressed by equation (9), with no specific relation required between 

U/U^ and • 

Comnarisons of the total shear stress as calculated from equation 
(11 ) with the measured distribution are included in figure 18, where the 
Ludwieg-Tillmann values of wall shear stress were used. The expermen- 
ta^oLl shear stress curves are within 10 percent (except at station 
3) o^ SLe predicted by equation (U), but at some stations the trends 
with y differ as the wall is approached. 

Figure 20 shows the total shear stress distribution for each sta- 
tion Such distributions are needed for accurate calculations 
bulent boundary layer development, where the integral relations of the 
e^SionHr^tion are employed (ref. 24). In the absence of such 




to 

rH 

to 
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distributions ; Granville (ref. 25) developed an expression for the inte- 
grated shear distribution through the boundary layer (this being the 
specific information needed to evaluate the relations of ref. 24) for a 
zero pressure gradient flow 'iItq. Evaluation of the sheajr distributions 
for adverse pressure gradients given in reference 6 led Granville to con- 
clude that the effect of pressure gradient was negligible^ namely 

t = ^ 0 - 

Comparison of the values of evaluated from the total stress 

curves of figure 20 , with the equation and the measurements of reference 
6 given hy Granville^ is shown in figure 21. A systematic variation in 
the meas\ored shear integrals of this report with the zero press\rre gra- 
dient values of is seen. However, the value at station 1, which is 

in a region of nearly zero pressure gradient, appears to agree well with 
the predictions. The variation indicates a definite effect of pressure 
gradient on t, which disagrees with the conclusion of Granville. It 
was noted, however, that the variation obtained herein was of no greater 
magnitude than the scatter appearing in the values of i|r evaluated by 
Granville from reference 6. 

Evaluation of mean V velocity . - Evalxiation of the momentum terms 
in the boimdary layer requires knowledge of the distribution of the 
component V (normal to the wall) of the mean velocity. The V com- 
ponent is usixally veiy small in the boundary layer, but it will appear 
in the momentum equation either as a derivative across the boundary 
layer or multiplied by some large quantity. Unlike the U component, 
the V component is not constant outside the boundary layer in a region 
of adverse pressure gradient. This result follows from the continuity 
equation together with the boimdary condition V=0 at y=0, which 
leads to 



The two-dimensional continuity equation in the form of equation 
(12) was used to detemine the values of V through the boundary layer. 
This form implies two-dimensional flow and thus neglects any lateral 
flow of the type suggested in reference 16. The mean V distributions 
obtained are shown as the solid curves in figure 22. 

A second method of evaluating the V velocities was from the flow 
deflection angles. The values of V calculated from the angles are in- 
cluded as the symbols in figure 22. Unfortunately, the two methods of 
evaluating V do not give comparable results except near the wall. 
Because of the small magnitude involved in the determination of V, 
either method may be subject to considerable inaccuracy; therefore it 
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was impossible to conclude which values best represent the actual quan- 
tities. The values of V obtained from the continuity equation were 
errployed for the analysis in the following sections, since a complete- 
set of measurements of the flow deflection angles was not obtained. 


Equation (lO) may also be used to evaluate the mean V component 
near the wall. The three methods have been compared in this region in 
figure 23, where the universal velocity profile was foimd to be valid. 

Stress tensor. - As noted in the section THEORETICAL BACKGROUND, 
one of the basic gro\;ps of terms encountered in a momentum study of vis- 
cous shear flows is the stresses. These terms form the components of a 
s-tness tensor (eq. which is composed of two additive tensors - 

pressure— viscous and turbulent terms . The terms of the turbulent stress 
tensor are of chief interest, since only limited information is avail- 
able on their distributions through the shear regions. 


The terms of the turbulent stress tensor for the flow investigated 


can be written 


— 



u^ 

uv 

0 

uv 


0 

0 

0 



(13) 


The measured values are plolled in figure 24. Xi may be seen 'bhai 

here again the normal stress pu^ is by far the largest term and also 
that the other three teirnis are of equal magnitude^ except near the "wall. 

As the pressure gradient increases^ pv^ and -puv approach the same 
value throughout the boimdary layer. This equality of -puv and pv^ 
was one of the basic ideas iinderlying Prandtl’s hypothesis of the mixing 
length theory (ref. 26) 


uv 



<3y 


(14) 


Although the mixing length theory is now generally considered to be in- 
adequately foimded, the assumption uv = v^ is nearly correct. 


4 . 


Also included in figure 24 are the stress terms due to the mean^ 
flow. Of these mean terms, it was found that the pressure contribution 
pA (see legend of fig. 24) to the normal stresses is by far the largest, 
but^his pressure contribution represents only the potential level of the 
system and, as concerns the actual forces, only the derivatives of the 
stress are of importance. 
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A second factor of interest in regard to the stress tensor is the 
direction along which the msiximum stress acts. For the turbulent stress 
tensor the angle 3 between the principal axis and the x-direction is 
given by the equation (ref. 27) 


tan 23 


2uv 



(15) 


The values calculated for 3 are shown in figure 25. Only the points 
near the wall can be evaluated accurately because eq\iatlon (l5) is 
extremely sensitive to the measured values in the outer edge of the 

boundary layer, since u^ -*■ v^ and uv -► 0. Disregarding the points 
for y > 1 inch (and excluding station 3), 3 was nearly constant for 
the different x-stations . As the wall was approached the angle decreased, 
so that the direction of principal stress was alined nearly parallel to 
the flow and the wall. 

Evaluation of terms appearing in equations of motion . - Information 
is now available from the present measurements to evaluate the terms 
appearing in the x- aind y-direction equations of motion, equations (2) 
and (3) . 


The terms of equations (2) euad (3) have been evaluated from the 
faired data of this report and are shown in figure 26, where the indi- 
vidual terms are multiplied by 9/U^ to give dimensionless values of 


n2u 

rate of change of momentum. The term v — - was fo\md in al 1 cases to 

be of the order of 1/1000 of v — ^ and therefore was not included in 

dy'^ 

the figures. Several theoretical investigators (refs. 14, 19, and 20) 

have speculated that shoxild be included in the equation, as it 

represents a contribution to the integral momentum equation; however, 
for the flows investigated, this term was found to be too small com- 
pared with the other terms of equation (2) to be plotted. 


Curves of the experimental difference between the left and right 
sides of equation (2) are also included in figure 26. This difference, 
which should be zero, is attributed largely to uncertainty in the term 


U which was quite sensitive to ciorve fairing. Only the balance at 
station 3 appears to show a large discrepancy, larger than the experi- 
mental error in ^ This discrepancy is thought to be due to the 

possible existence of two types of equilibrium profile during the meas- 
urements of shearing stress and mean velocity. 
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The resiilts of the evaluation of the terms of equation (s) appeared 
at first to be rather unexpected. The usual assimiption made in boundary- 
layer studies is that this equation may be neglected completely, and so 

it might be, except for the turbulent transfer velocity v^. All the 
mean flow terms of equation (3) are small and might well be neglected, 

but the term was found to be as large, if not larger, than any of 

the terms appearing in the x-direction equation of motion. At firs-t there 
appeared to be a discrepancy in the terms evaluated for ^uation (3), as 


Sv2 

no other term appeared capable of canceling the large term. A re- 

examination revealed that a small \mcertainty in the measured values of 

-j 

static pressure would result in a relatively large value of p, If 
now it is assumed that the only terms of appreciable magnitude in equa- 
tion (3) are ^ and i equation (3) may be integrated to give 


pv^ 


= PW - P 


(16) 


since pv^ = 0 and p = P^ at y - 0. 

The magnitude of change in static pressxxre can be estimated from 

equation (l6). A check of the data shows the maximum value of pv^ 
measured was approximately 0.008 pound per square foot (fig- H)^ which 
results in an absolute static -pressure variation of roughly 0.0004 per- 
cent. During the course of meas-urements , some static-press-ure surveys 
were made, with indications being that any variation in static pressure 
was of the same order of magnitude as the uncertainty in the measure- 
ment. The variation in static pressure as given by equation (16) was 
found to fall within the region of uncertainty of the measirrements . 


For fully developed incompressible turbulent channel flow the y- 
momentum equation (3) is known to reduce exactly to equation (16) (ref. 
22) and for the present case of turbulent boundary layer flow, the same 
reduction appears to be a good approximation. This last conclusion is 
also reached in reference 28. Comparison of the x-derivative of equa- 
tion (16) with equation (2) shows that a^ coupling between the two 

^2 

equations must reside in the vulue of experimentally, this coupling 


was 


found to be negligible the x- and y-momentum 

equations (2) and (3) are concluded to be independent for the flow in- 
vestigated despite the presence of comparably large terms in both. 
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Evaluation of energy terms available from measurements . - The ainal- 
ysis of this report has been directed primarily toward the evaluation of 
the terms important in the study of the momentum distribution of the 
flow; however^ a brief analysis of the distribution of the kinetic 
energy of the turbulent motion has been made. 


The distributions of the local turbulent kinetic energy per imit 

mass 


E 


u^ 


+ v^ + w^ 


local 


(17) 


are presented in figure 27 . Of more interest is the average energy con- 
tained in the turb\ilent motion, which is written nondimens ionally as 


E. 


'total 





Elocal 


(18) 


The resulting values of E^otal shown plotted in figure 28. Start- 

ing with the turbulent energy equation (first moment of the momentum 
equation)^ the author of reference 29 was able to arrive at an integral 
equation from which the value of might be estimated. The anal- 

ysis of reference 29 was made for fully developed flow in pipes and 
channels^ with the resulting evaluation giving the equation 


which gave reasonable correlation with the fully developed channel flow 
of » reference 3. Equation (l9) was fitted at the first x-station and the 
resulting curve is included in figure 28. Thus it would appear that the 
predictions of reference 29 developed for the one -dimensional fully de- 
veloped turbulent flows are not incompatible with the turbulent energy 
in the boundary layer flows with arbitrary pressure gradients 
investigated. 


The rate of removal of kinetic energy per unit volume from the mean 
flow by the turbulent stresses is (ref. 6) 


"T dU SV 

u^ — + V ^ + 


dx 


Sy 


uv 


/SV 

\bx ^ Sy/ 


The distribution of this terni through the boundary layers investigated 
is presented in figure 29. This energy derived from the mean flow goes 

— SU 

directly into the production of turbulent motion. The term uv 

which is equivalent to t was by far the largest term over most of 
the boundary layer. 
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CONCLUSIONS 


Analysis of measurements in turbulent boundary layers with adverse 
pressure gradients has led to the following conclusions: ' 

Measurements show that the turbulent velocities and are 

nearly equal over most of the boundary layer, except near the wall, and 
that they are roughly one -half of the longitudinal velocity fluctuation 

The shape of the curve was foiind to change with increasing 

pressure gradient. Measurements indicate that two possible distributions 
coilLd exist, depending on flow conditions. A change with pressure gra- 
dient was noted for the mean velocity profiles. 


At all pressure gradients good agreement was obtained for the values 
of skin friction determined from the Ludwieg-Tilljiiann relation and those 
measured with the heat transfer - skin friction instrument. At the low 
pressiire gradients the agreement extended, in addition, to the Falkner 
flat -plate equation. 


Use of the measured distributions of turbulent shear stress -puv 

and viscous shear stress p to predict the wall shearing stress 

agrees roughly with the values obtained from heat -transfer measurements 
or predicted by the Ludwieg-Tillmann equation. 


Measurements indicate that it may be possible to assume a sublayer 

within which the total shear stress obeys the relation ^ where 

T is the total shear stress eind py is the static pressure at the 
wall. 


The universal relation 

mean velocity, is the shear stress velocity, and v is the kinematic 
viscosity) that has been observed by many workers was found to hold near 
the wall, and the total shear stress distributions predicted by use of 
this universal relation were found to agree qualitatively with the dis- 
tributions obtained from the measirrements . 


against 




(where U is the local 


The variation of the integrated shear distribution through the tur- 
bulent boundary layer suggested by Granville for zero pressure gradient 
does not appear to hold for the adverse press\ore gradients investigated 
herein, although for the one station where the pressure gradient was 
nearly zero the agreement was good. 
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Evaluation of the terms of the turbulent stress tensor indicates 

that the normal stresses pv2 and pw2 and the shear stress -puv are 

of the same order of magnitude^ with pv^ and -puv being nearly equal 
over the whole region for the largest pressure gradient. 


The angle between the principal axis of the turbulent stress tensor 
and the x-direction was found to be nearly constant with x-distance and 
independent of press\rre gradient. The angle also shows that the direc- 
tion of principal stress near the wall tends to be alined nearly parallel 
to the wall. 


Evaluation of the terms in the x-direction equation of motion indi- 
cates the equation may be written to sufficient accuracy as 


U 


Su 


+ V 


Su 


Suv 


+ V 




where V is the mean velocity component normal to the boundary^ p is 
the air density^ and p is the static pressure. All other terms are 
small for the flows investigated. 

Evaluation of the terms in the y-direction equation of motion indi- 
Sv^ 

cates the term -g- — to be as large as or larger than any of the terms 

in the x-direction equation. A sirfficient approximation to the y- 
direction equation of motion appears to be 


_ 1 Sp 

by ~ ~ Q by 

or 

pv^ = - P 

for the flows investigated. The x- and y-equations of motion are inde- 

^v2 1 rVrj 

pendent if «i which is true for the cases investigated, 

ox p ox 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland^ Ohio, August 10, 1954 
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APPENDIX A 
SYMBOLS 


The following symbols are found in this report: 

A constant, depending on physical quantities of hot wire 

B constant, depending on physical quantities of hot wire 

local wall shear stress coefficient, 

T 


1/2 pU^2 


E 


local 


E. 


total 


As 


H 


I 

P 

P, 


-W 


R 

R 


shearing stress coefficient, — ■ ^ 

1/2 

local turbiilent kinetic energy per unit mass 
average energy contained in turbulent motion 
voltage output of hot wire 

fluctuation voltage output of hot wire due to mass -flow 
f luctuat ion 

mean- velocity -profile form parameter, 6*/9 

current flowing through hot wire 

mixing length 

static pressure 

reference static pressure 

static pressure at wall 

operating resistance of hot wire 

resistance of hot wire at flow temperature 


a 

Re 


Reynolds number based on a convenient length; for this 

pUqS* 
report, 


Re Reynolds number based on momentum thickness. 


pU.0 


3173 


hot-wire sensitivity 
stress tensor 
turbulent stress tensor 

local mean velocity in x- and y-directions 

mean velocity of free stream where viscous effects are 
neglected 

shear stress velocity^ 

instantaneous resultant mean plus turbulent velocities 

instantaneous turbulent velocities in x-^ y-, and z-directions 

mean squares of turbulent velocity fluctuations in x-, y-^ 
and z-directions 

Reynolds turbulent shear stress in xy-plane 

Reynolds turbulent shear stress in xz -plane 

Reynolds turbulent shear stress in yz -plane 

direction parallel to test wall and in general direction 
of mean flow 

direction normal to test wall and approximately normal to 
mean flow 

direction normal to rneatn flow and paxallel to test wall 

angle in xy-plane resultant mean flow makes with wall 

angle between principal axis of turbulent stress tensor and 
x-direction 

boundary layer thickness 

displacement thickness^ representative of mass -flow deficiency, 



Krone cker delta symbol 
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P- 

V 

P 

2 

T 

T 


w 


3{X 


yy 


'ZZ 


xy 


<P 


momentum thickness, representative of mean -flow momentum 
deficiency, / if-ll - — jdy; corrected for turbulence^ 


/ 4 (‘ ■ -)■ 


rP 


u 

1 

*^xx 


Up 

L ""A 

pU^ 

pu/. 


dy 

, I I I 1 _ % r. / 

/O 

viscosity of air 

kinematic viscosity of air, n/p 

density of air 

instantaneous angle between x -direct ion and instantaneous 
velocity vector 

local shear stress 

local shear stress at wall in x-direction 
normal stress terms, defined in eq. (4b) 


^ 4 . 4 . — 

shear stress term, 'Sxj " 


angle between x-direction and hot wires 

^6/6* 


0 


integrated shear distribution, ^ 2c^ \^j 

integrated shear distribution for zero pressxrre gradient 
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APPENDIX B 


THEORY OF CONSTANT -TEMPERATURE HOT-WIRE ANEMOMETRY MEASUREMENTS 

The constant -temperature hot-wire anemometer is a system designed 
to eliminate the problem of thermal inertia enco\mtered in the more 
common constant-current hot-wire anemometer systems. The use of the 
d-c feedback amplifier system employed in the measurements of this re- 
port to correct for any change in the wire operating temperat\jre was 
suggested in reference 30. 

The constant -temperatirre system employs the d-c amplifier as a 
servomechanism to sense the unbalance of a Wheatstone bridge, of which 
the wire is one arm, and by continuously varying the feedback voltage 
the bridge is kept in balance. If the wire is the only variable resis- 
tance in the bridge , the feedback voltage will serve to keep the wire 
resistance and likewise the wire temperature constant; thus the problem 
of thermal lag of the wire is eliminated since the change in wire tem- 
peratiire is negligible. A measure of the fluctuating feedback voltage 
is an indication of the fluctuating mass flow over the wire. 

The refinements of the constant -temperat\ire system of this report 
over that appearing in reference 12 are the shock mounting of the ampli- 
fiers and the replacing of the electronic tube rectifiers of the ampli- 
fier power supply by selenium rectifiers . Both changes were made to 
eliminate as much amplifier backgroimd noise as possible. The selenium 
rectifiers also helped to relieve a problem in amplifier cooling, which 
in turn helped to eliminate noise . 

The output of the system was read on the average -square computers 
described and evaluated in reference 12. The sum and difference circuit 
was also the one described in reference 12 in connection with the 
double -correlation instrument. The sum circuit gave accurate results, 
but as will always be the case the difference circuit was affected by 
the circuit noise level. If the same sine wave of 1 volt was fed into 
the two inputs and subtracted through the circuit, a difference mil l of 
approximately 5 percent of the input signal was obtained. The actual 
error encountered because of the difference null was negligible except 
near or in the free -stream t\irbulence, and fortunately here the contri- 
bution of the sum and difference tenns (see appendix C) to the total 
velocity fluctuation term was negligible. 

The use of a d-c amplifier requires that the bridge of the constant - 
temperature hot-wire anemometer system operate at a slight imbalance; 
otherwise the amplifier would oscillate. The unbalance of the bridge 
was maintained constant at 3 microamps by use of a microammeter. This 
unbalance results in a slight systematic error in the fluctuating feed- 
back voltage output . However, the data were not corrected as the ran- 
dom scatter in the measurements was found to overshadow the smaller 
systematic error. 
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The empirical relation between heat loss and mass flow for small 
cylinders electrically heated may be written 

.2^ 

= A + (Bl) 


where i is the operating current of the hot wire, R is the operating 
resistance of the hot wire, R^ is the wire resistance at some refer- 
ence temperature (usually that of the surrounding air stream), pU is 
the mass flow past the wire, and A and B are empirical constants 
which must be determined for each wire. Equation (Bl) has been referred 
to as King's equation (ref. 3l) since he first arrived at the form from 
a theoretical analysis of the heat loss from small circular cylinders in 
potential flows. However, the equation is justified from empirical re- 
sults, since King’s analysis did not take into accoiint viscous effects 
on the heat transfer. 

Replacing the curve of i against U as predicted by equation 
(Bl) by its tangent at any point makes possible the prediction of the 
instantaneous heat loss from small wires. Schubauer and Klebanoff 
(ref. 32) indicate that the heat loss is a function of the component of 
mass flow normal to the wirej thus equation (Bl) can be generally 
written 


R 




- R 


a 


= A + B-^pVj^ sin 2 


(B2) 


where V. is the total instantaneous velocity and Z is the instan- 
taneous iigle between and the wire. Differentiating equation (B2) 

with respect to the velocity vector, where the wire resistance is held 
constant and only the wire current varies, gives (incompressible flow) 


2iR di ^ B P 

R - Rq, 2 sin Z 


(sin Z dV^ + cos Z 


dZ) 


(B3) 


'pjixs 0 qu 3 ,tion will be true for small fluctuations, where only the linear 
terms of the fluctuations are retained. In order to rearrange equation 
(B3) into a form where hot-wire measurements can be analyzed, the fol- 
lowing approximations are made: (l) For small changes in Ae, the fluc- 

tuating voltage, let R di = de = Ae. (2) Assume the instantaneous 
velocity is made up of a mean component U which makes an angle (p 
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with the wire (see sketch l) plus the turbulent fluctuation components 
u and V. For this derivation is taken as coplanar with U and 

the wire. With the restriction that u and v are very sma i i compared 
with it follows from sketch 1 that 

(a) dV^ “ u 

(b) dZ — 

(c) sin Z - U sin <P 

(d) Z == (p 

Equation (B3) can now be written 

u sin (p ± V cos (P = - Ae (B4) 

CL 


Equation (B4) is the general equation relating the fluctuating voltage 
due to the d-c amplifier feedback to the turbulent velocity fluc- 

tuations occurring across the hot wire. This derivation parallels that 
given in reference 32 for the constant-current system. For the simplest 
case^ that of a wire normal to the mean velocity equation (B4) re- 
duces to 



The temns on t he right side of equation (B5) are determined at each 

point where is measured. The current is determined by measuring 

the mean bridge current and then using a balanced Wheatstone bridge re- 
lation to calculate the current through the wire. The mass flow VpU 
and B are obtained fro)r a calibration curve for the wire. The wire 
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is placed in an air stream where the mass flow can he varied and deter- 


mined by other means. A plot of the mean heat loss 


i% 

(R - Ra) 


against 


the square root of the mass flow ■V^ is obtained. From this plot B 
is obtained as the slope of the curve ^ and can be determined from 


Jc\ 

the value of measured at each point. The resistance R can be 

R - K 

(X \ 

obtained from Wheatstone bridge relations, and the value of is 

usually measured directly by use of a potentiometer. The fluctuating 
voltage was measured in terms of its root mean sq\iare value by means of 
an electronic squaring circuit. 


The evaluation of hot-wire data in several cases required correc- 
tions for changes in heat loss calibration in order that consistent 
results might be obtained. The initial procedure was to calibrate the 
wires in the free stream of tne channel after a profile through the 
boundary layer was completed. This type of calibration resulted many 
times in discontinuous jumps observed in the calc\ilated quantities. 

(The discontinuous jumps were apparently caused by dirt particles strik- 
ing the wires.) It was found that consistent results could be obtained 
if the velocity profiles measured in the boundary layer with a total- 
pressure probe were also employed. From the measured mean veloci^ pro- 
file for a given station and free-stream velocity, a value of 
each point in the boundary layer was determined. If now the heat loss 
of the wires at each point in the boundary layer is plotted against the 
predetermined -y/ pU, a master calibration of the wires at each point is 
obtained. In general, it was found that changes in calibration resulted 
only in shifting the calibration curve, with the slope remaining nearly 
constant . 
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APPENDIX C 


REDUCTION OF x-WDRE DATA 


The general equation relating the voltage fluctuation output of the 
anemometer to the turbulent fluctuations was derived, in appendix B. 


u sin (P ± V cos (p = Ml/pU ^ 

pB (R-R^) 


(B4) 


The term sin (p 


pB(R-^ ) — sensitivity factor S of the wire^ since 

it depends only on the mean quantities; thus the equation can he written 


u sin <P + V cos <P = S As 


(Cl) 


If two wires are placed in the flow so that the mean velocity U makes 
an angle <p with each, as shown in sketch 2, it may be seen that 



Sketch 2 

equation (Cl) gives the following eqimtions for the two wires ; 


u 


sin cp + V cos (p = 


and 


sin^ <p + 2uv sii 

,2 r,-;-r»2 /ft _ o7777 


Sj_Aei 

(C2) 

^2^2 

(C3) 

1 are 


' cos^ cp = 

(C4) 

cos^ cp = Sg^Aeg^ 

(C5) 
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Taking the difference of equations (C4) and (C5) yields 


uv = 


4 sin (p cos cp 


(C6) 


which is the relation used to determine the turbulent shearing stress 
puv. 


The sum of equations (C4) and (C5) gives 


u*^ sin^ <p + cos ip = 2 


(C7) 


In order to obtain another equation between u^ and it is neces- 

sary to obtain electrically the instantaneous sums and differences of 
the voltages of equations (C2) and (C3) . The mean square of the instan- 
taneous sum may be written 

■- ■ '2 2 2 

sin^cp(S2+Si)^ + 2uv sin (p cos <p (S2-Sx) (S2+S1) + v cos <p(S2-Si) 

= S^^^Sg^CAe^ + Ae2)^ (C8) 

The mean square of the instantaneous difference may be written 

u^ sin^<p(S2-S-L)^ + 2uv sin «P cos <p (82-33^) (82+82^ ) + v^ cos^cp(82+Sx)^ 

= - A=2)^ (C9) 

The difference of equations (C8) and (C9) yields 

^ sin2(P - ^ cos2(p= (CIO) 

where 

e^2 = (Ae^ + Acg)^ 

t} (Aei - Ae2)^ 
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Solving equations (C7) and (CIO) for and gives 



(Cll) 


and 



2 _ e 2 


(C12) 


It should be noted that v can be replaced by depending on the wire 
orientat ion. 

Figure 30 shows a block diagram of the instioimentation hookups used 
to obtain the turbulence measurements of this report. 
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Figure 1. - Schematic diagram of 6- by 60- inch subsonic boundary- layer channel 
(All dimensions are in inches.) 
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Figure 2. - Static -pressure distribution along tunnel wall. Reference static pressure p^, 
6.600 pounds per square foot (p^ = atmospheric pressure - tunnel static pressure at 
X = O) (atmospheric barometric pressure, 29.21 in. Hg. at 81.5^ F). 
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Figure 3. - Pressure and hot-wire boundary- layer probes. 
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Figure 4. - Voltage control circuits for pressure transducers. 
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Figure 5. - Total-pressure variation ot)served through boundary layer at 
station 1. Each division on time scale equals approximately 10 seconds j 
smallest division on pressure chart is approximately 0.120 pound per 
square foot. 
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(a) Station 1. feet per second. 

Figure 6. - Mean velocity distrilDution through boundary layer. oj 
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(■fa) Station 2. * 52 feet per second. 

Figure 6. - Continued. Mean velocity distribution through boundary layer. 
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(c) Station 3 . • 43 feet per second. 

Figure 6 . - Continued. Mean velocity distrilDution through boundary layer. 
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(d) Station 4. U]_ » 43 feet per second. 

Figure 6. - Concluded. Mean velocity distribution through boundary layer. 
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Figure 7- - Boundary- layer flow angle measuring probe. 
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Figure 8. - Mean flow deflection angle through "boundary layer. 
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Figure 9. - Variation of local skin-friction coefficient along test wall. 


cn 


NACA TN 3264 


Longitudinal turbulent intensity, 


46 


NACA TN 3264 


.10 


.06 


.04 


^ -02 














I 









Probe 

□ 

A -/v^ 
o 

\ 










Sal 

) 









<r 

o< 
















i D 

2 

□ 












4t 

} 













r> 

□ 





■1 







i 

* □ 

o 

1*1, A 



3 

G 

L 

6* 

1 








u A 

[ 

I □ 

n 

□ 


.25 


.50 .75 1.00 1.25 1.50 

Distance from wall, y, in. 

(a) Station 1. 


.10 


.08 


.06 


.04 


.02 












' 


S— 

o 

























°0 












— ^ 

( 

) 

o 












( 

> 








O 




O 

< 

) 







O 

B 





( 

> 












{ 

) 

< 

> . 

> 1 

) 


6 

1 

1 5* 

J 












.4 .8 1.2 1.6 Vi.u 

Distance from wall, y, in. 

(b) Station 2. 

Figure 10. - Distribution of longitudinal turbulent intensity. 
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Figure 10. - Concluded. Distribution of longitudinal turbulent intensity. 
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Figure 12. - Distribution of 
through boundary layer. 
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(a) Station 1. 



(h) Station 2. 

Figure 13. - Turbulent shear stress distribution through boundary layer. 
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(d) Station 4. 


Figure 13. - Concluded. Turbulent shear stress distribution through 
boundary layer. 
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Figure 14. - Determination of uv from outputs of two hot wires at 45 to mean flow. Station 2. 
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Figure 15. - Variation of mean-velocity-profile parameters along test wall. 
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Figure 16. - Variation of mean-velocity-profile form parameter and Reynolds number along test wall. 
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(a) Station 1. 

Figure 17. - Distritiution of mean and fluctuating velocities through boundary 
layer. 6 = 0.136 inch) 0 = 0.103 inch) H = 1.32) Reg = 2740. 
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(*b) Station 2. 

Figure 17. - Continued. Distribution of mean and fluctuating velocities through 
boundary layer. 6 = 0.225 inch; 6 = 0.169 inch; H = 1.33; Re^ = 4377. 
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(c) Station 3. 

Figure 17. - Continued. Distribution of mean and fluctuating velocities through 
boundary layer. 5 = 0.369 inch; 0 = 0.257 inch; H * 1.43; Re^ = 6146. 
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(d) Station 4. 

Figure 17. - Concluded. Distribution of mean and fluctuating velocities through boundary 
layer. 5 =0.652 inch; Q = 0.393 inch; H = 1.65; Re^ = 7800. 
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(a) Station 1. 

Figure 18. - Shear stress distribution near wall. 
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Figure 18. - Continued. Shear stress distribution near wall. 
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(c) Station 3 . 

Figure 18. - Continued. Shear stress distribution near wall. 
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(d) 'Station 4. 


Figure 18. - Concluded. Shear stress distri'bution near wall. 
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Figure 19. - Semilogarithmic plot of mean velocity distribution near wall. 
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Figure 20. - Total shear stress distribution through boundary layer. 
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Figure 22. - Variation across boundary layer of mean velocity component normal to boundary. 
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Figure 23. - Predicted values of mean V velocity near wall. 


3173 


MCA TN 3264 


67 


to 

t-- 


to 


M 

0 
(d 

rQ 

CD 

1 

PR 

o 



(b) Station 1; mean motion. ^ = 178,700 
Figure 24. - Distribution of stress terms through boundary layer. 
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Figure "4. - Continued. Distribution of stress terms through boundary layer. 
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(e) Station 3; tur^bulent motion. 



(f) Station 3; mean motion. ^ = 308,900 
Figure 24. - Continued. Distribution of stress terms through boundary layer. 
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Figure 24. - Concluded. Distribution of stress terms through boundary layer. 
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Figure 25. - Variation through boundary layer of angle between principal axis of turbulent 
stress tensor and x-axis . 
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(b) Station 1; y-dlrectlon. 

Figure 26. - Distribution of momentum terms through boundary 
layer. 
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' Figure 26. - Continued. Distribution of momentum terms through 
boundary layer. 
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(f) Station 3; y-direction. 

Figure 26. - Continued. Distribution of momentum terms through boundary 
layer . 
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(h) Station 4; y-direction. 

Figure 26. - Concluded. Distribution of momentum terms through boundary layer. 
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Figure 27 . - Turbulent energy distribution through boundary layer . 
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Figure 28. - Variation of "bulk turbulent energy along test wall compared with prediction of 
Corrsin for a fully developed channel flow. 
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Figure 29. - Rate at which kinetic energy of mean flow is transformed into 
turbulent kinetic energy. 


3173 


3173 


MCA TN 3264 


79 




I 



Constant- temperature 
hot-wire anemometer 
amplifier 


(a) 




Measurement , 





Constant-temperature 
hot-wire anemometer 
amplifiers 

(h) X-wire measurments. 


CD- 3753 


Figure 30. - Block diagram of instrumentation hookup. All connecting cables 
are of the coaxial type. 
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